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Abstract: Black holes are fascinating objects. As a class of solutions to the Einstein 
equations they have been studied a great deal, yielding a wealth of theoretical results. 
But do they really exist? What do astronomers really mean when they claim to have 
observational evidence of their existence? To answer these questions, I will focus on a 
particular range of black-hole masses, approximately from 10 2 M Q to 10 4 M Q . Black 
holes of this size are named Intermediate Mass Black Holes (IMBHs) and their exis- 
tence is still heavily disputed, so they will be perfect for illustrating the observational 
challenges faced by a black hole hunter. 



Sazetak: Crne rupe su fascinantan objekt. One su, kao klasa rjesenja Einsteinovih jed- 
nadzbi, bile predmetom mnogih istrazivanja koja su dovela do obilja teorijskih rezul- 
tata. No postoje li one zaista? Sto astronomi tocno misle kada turde da imaju opser- 
vacijske dokaze da one postoje? Da bih odgovorio na ova pitanja usredotocit cu se na 
odredeni raspon mase crnih rupa, onaj izmedu 10 2 M e i 10 4 M Q . Crne rupe ovih di- 
menzija nazivaju se crne rupe srednjih masa (IMBHs) i njihovo je postojanje jos uvijek 
nasiroko osporavano, tako da ce one biti savrsene za ilustriranje opazackih izazova s 
kojima se suocava lovac na crne rupe. 



Sommario: I buchi neri sono oggetti affascinanti. Come classe di soluzioni delle 
equazioni di Einstein, sono stati studiati approfonditamente, producendo una messe di 
risultati teorici. Ma esistono davvero? Che cosa intende dire un astronomo quando 
afferma di avere delle prove osservative della loro esistenza? Per rispondere a queste 
domande, mi concentrero sull'intervallo di masse approssimativamente compreso tra 
10 2 M Q e 10 4 M e . I buchi neri di queste dimensioni vengono chiamati buchi neri di 
massa intermedia (IMBHs) e a tutt'oggi si discute ancora se esistano o meno. Saranno 
quindi perfetti per illustrare le difficolta osservative incontrate da un cacciatore di buchi 
neri. 



1 Introduction 



Intermediate Mass Black Holes (IMBHs) have masses in the range 10 2 - 10 4 M . 
Stellar-mass black holes are the product of normal stellar evolution. Tipic ally they 



weigh several solar masses but some have been found heavier than 10 Mn (seelCreighton & Price 



20081) . To the best of my knowledge the record to date is 15 M Q (Bulik 2007). The so- 
called supermassive black holes, instead, weigh billions of solar masses and have been 
found in the centers of galaxy bulges, and they are thought to power much of the high- 



energy activity going on there, (see iKormendv & Richstond[l995t iFerrarese & Ford 



2005, for reviews). IMBHs live in-between these two extremes and are actually a 
promising bridge between them, as they are speculated to constitute the building blocks 
of supermassive black holes. 

Before plunging into the mysteries of IMBHs we shall deal for a while with stellar 
mass and supermassive BHs. 

1.1 Stellar-mass black holes 

The formation mechanism of stellar-mass black holes is reasonably understood, as it 
is based on the time-tested theory of stellar evolution. Black holes are basically what 
is left after stellar death. How a star dies is dependent on its mass at the time the 
nuclear reactions capable of sustaining it against the pull of gravity stop. During the 
giant-branch phase of its evolution, a star typically looses mass in variable amounts, so 
its final mass is not precisely constrained. Anyway, when nuclear reactions stop, the 
matter constituting the star collapses, unless other supporting mechanisms intervene. 
Electron degeneracy pressure is one of these mechanisms, and if the mass of the col- 
lapsing star is small enough it will win over gravity and the star will become a white 
dwarf. If electron degeneracy pressure yields to gravity, the matter is converted into 
neutrons and the collapse continues. Neutron degeneracy pressure is the next barrier: 
if the mass is not high enough, the collapse is stopped by it and a neutron star forms. 
If the mass is high enough, no known mechanism can stop the collapse and the star 
is expected to shrink until it is contained within its own Schwarzschild radius, thus 
becoming a black-hole. Most evidence about the existence of stellar-mass black holes 
comes from X-ray binaries. In these objects, an accretion disk forms and becomes hot 
enough to emit in the X-ray band (see lMcClintock & Remillardll2003l) . X-ray emission 



plays a crucial role in the context of IMBHs as well, and will be covered in more detail 
later. 



1.2 Supermassive black holes 

Supermassive black holes were first postulated as engines ca pable of sustaining the 
high-power emission of quasars (e.g. IFerrarese & Fordl 120051) . In the local universe 
they have been observed in the bulges of several galaxies either thanks to the accreting 
hot gas or indirectly by their effects on stellar kinematics. Also in the bulge of our own 
Galaxy, the Milky Way, a supermassive black hole was found. It is the supermassive 
black-hole with the best measurement of mass to date, as the orbi ts of resolved sta rs 
were used to estimate the black hole mass with w 10% uncertainty ( Ghez et al] l2003). 
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1.3 Why focus on IMBHs? 

We have already stated that IMBHs naturally bridge a gap between two observational 
regimes where relatively firm detections exist. There are further reasons that make this 
class of black holes interesting. 

First of all, IMBHs are likely good candidates as gravitational radiation sources. If they 
are to be found in a large fraction of massive star clusters (see below why it is thought 
that star clusters might contain IMBHs), they are likely hundreds of objects in the Local 
Group alone. In the dense environment of a star cluster, close encounters with stellar- 
mass black holes and other massive objects would be quite fre quent, leading to a go od 
potential for gravitational wave emission. Estimates for LISA (IFregeau et alj|2006l) or 
a similar space-borne interferometric detector of gravitational radiation predict tens of 
IMBH in-spiral events per year. Somewhat lower nu mbers are ex pected for ground 
based interferometers such as VIRGO or LIGO (see lMillerj|2002l) . which work in a 



different frequency range and would preferentially reveal black-hole merging events. 
In a typical galaxy, thousands of star clusters which could potentially contain an IMBH 
are present. Dynamical friction due to gravitational encounters with stars could lead to 
an infall of star clusters (with the associated IMBHs) to the galactic bulge. Actually, so- 
called nuclear star clusters have recently been observed in the bulges of several galax- 
ies. If these where to contain IMBHs they could provide a way to see d super-massive 



black holes, suggesting an explanation for the formation of the latter dEbisuzaki et al 
200ltlVolonteri et alJl2003l:lTanaka & Haimanll2009l) 



More generally, IMBHs would have a strong effect on star-cluster dynamics, as they 
can play a crucial role in the gravothermal collapse of its core (core collapse), basically 
halting it. Core collapse is one of the processes resulting from the collisional evolution 
of clusters of stars. Stars in a star cluster evaporate as they diffuse in phase space, at- 
taining velocities that exceed the cluster escape velocity. This way, energy is removed 
from the cluster, which contracts its core and moves its most massive stars to the center, 
in order to balance this effect (mass segregation). The three processes of evaporation, 
core collapse and mass segregation are all linked together in the dynamical evolution 
of a cluster. The amount of mass segregation attained b y stars in a relax ed star cluster 



is also heavily dependent on the pres ence of an IM BH (Gil l et alj|2008l) . which, inci- 
dentally provides a detection method (Pasqu ato et al.ll2009b . 



The density and velocity dispersion cusps ( Bahcall & Wolf 19761 : Novola & Gebhardl 



2006, 2007) that have been observed in many star clusters could also be explained by 
the presence of an IMBH. 

Finally, IMBHs would also naturally explain the so-called Ultra-Luminous X-ray sources 
(ULXs), X-ray emitting object s which appear t o radiate above the Eddington limit for 
a stellar-mass black hole (e.g. Makishi ma et alHoOOl) . 



2 IMBH formation 

Several scenarios for forming black holes in the 10 2 - 10 4 M Q mass range have been 
proposed. 

• Population III stars, i.e. the first stellar structures to form from primordial- 
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composition gas after the Big Bang, are predicted to be very massive. Their 
fast evo lution can end with a black-hole remnant o f several hundreds of solar 
masses dMadau & ReeshooUlSchneider et alll2002l) . 

Simulations of dense clusters with young stars sh ow that they may merge in a 



runaway cascade, producing a very massive star (Porteg ies Zwart et 



Such a star may produce a large black hole when it dies, 

Lower-mass black holes may merge and p roduce IMBHs in suffic iently high 
density environments. Dense star clusters (Mi ller & Hamiltonll2002h may thus 



produce IMBHs, if the accreting object is not kicked out by gravitational wave 
recoil, 

Primordial black holes in the intermediate-mass range may form in the early 
evolution of the univers e, according to certain variants of inflationary cosmology 
dKawaguchi et al.ll2008h 



Two out of four mechanisms involve dense environments, so the search for IMBHs 
is tied to the physics of star clusters in a natural way. Two or more formation routes 
can be active together, leading to a composite IMBH population. 

2.1 Formation from population III stars 
2.1.1 Cooling and fragmentation of primordial gas 

Population III stars are extremely metal poor, because they are formed from primor- 
dial composition gas. It is the physics of this gas, in particular the cooling mech- 
anisms, which determines how discrete clumps of gas form and collapse to become 
population III protostars. How the gas cools sets the final mass of the clumps pro- 
duced, eventually set ti ng the mass of population II I stars. Studies on this topic dHovld 
1953HMestellll965albl: Low & Lvnden-Bell 1 19761 ) show that cooling is inefficient in 



primordial-composition gas, which implies that population III stars are characterized 
by much higher masses than those of stars forming in the local universe. With stellar 
masses as high as hundreds of solar masses, IMBHs can easily be produced as rem- 
nants. 

A gas chunk of mass M undergoes gravitational collapse if the typical thermal 
energy kT of a particle is smaller than its binding energy: 

GMm 

kT < (1) 

R 

where R is the size of the chunk and m the particle's mass. Introducing the density 
p « M/R 3 we obtain 

M> Mjx {kTIGmf l2 p- lil (2) 

where Mj is the Jeans mass of the gas for a given density p and temperature T. A 
gas cloud of mass M » Mj will split into smaller fragments of mass ~ Mj, which 
collapse separately. As each chunk contracts, its density and temperature change, and 
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so does Mj(T,p). Increases in Mj lead to no further fragmentation, while a decreasing 
Mj makes the collapse of smaller chunks possible. 

If the chunk of gas is unable to efficiently cool, the collapse leads to a temperature 
increase. The extreme case is adiabatic collapse, where 

T oc p y- 1 (3) 

where y = c p /c v = 5/3 in the case of a purely monoatomic gas. By substituting Eq.[3] 
into Eq. [2] we obtain 

ocp I/2 (4) 

where the last relation holds for y = 5/3. Density increases during collapse and Mj 
increases with its square root, therefore no further fragmentation takes place. 

The opposite case is that of isothermal collapse, where any heat produced by col- 
lapse is immediately dissipated. Equation|2]then yields 

Mj^p- l!1 (5) 

so the Jeans mass is decreasing as collapse proceeds, and further fragmentation is al- 
ways possible. The more efficient the cooling mechanisms, the smaller the mass of the 
clumps of gas produced at the end of fragmentation. 



2.1.2 Gas cooling mechanisms 

Understanding the physics of gas cooling is crucial to determine the typical mass of the 
protostars that will form. A collapsing cloud of gas relies mainly on radiative cooling to 
eliminate the heat produced by the collapse. Radiative cooling depends on the radiative 
de-excitation of atoms and molecules excited by collisions. 

Solar-composition gas contains metals, which are extremely efficient coolants. To 
name a few atomfl C II, Ne II, S i II, N II, O II, O III, have transitions which span 
the energy range 2.5 x 1CT 3 eV - 5.35 eV. This corresponds temperatures from « 30 
K to « 6.2 x 10 4 K, which insures the gas is able to efficiently cool over a wide range 
of te mperatures. This entail s a decreasing Jeans mass until about 1M clumps form 
(e.g. lBromm & Larsonll2004 . Primordial gas contains virtually no metals, and thus is 



must cool through less efficient channels. Also dust, which plays a role in the cooling 
of solar-composition gas, cannot be formed in the absence of metals. The current un- 
derstanding is th at molecular hydrogen, H2, is the main coolant in primordial gas (e.g. 
Brommetal. 2002). This comparatively inefficient cool i ng leads to protostar masses of 



1 00 M ( Bro mm et all2002t iBromm & Larsonll2004 lGloverll2005t lNormanll2 008 ) . 



2.1.3 Evolution of primordial-composition stars 

As shown, population III stars are predicted to be very massive based on the physics 
of primordial gas cooling. Do they loose lots of mass during their evolution, or can 
we trust them to keep it until death, allowing the birth of an IMBH? The massive O- 

1 In the following we make use of the spectroscopic notation, where A I indicates the neutral atom of 
element A, A II the singly ionized atom, A III the doubly ionized, and so on. 
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and ff-type stars in the l ocal universe lose mass through radiatively driven winds (e.g. 
Lucy & Solomoiilll970t ICastor et alJll975t ICassinellilll979t lAbbottlll980i) . The force 
per unit mass / exerted on the gas in the outskirts of such a star depends on the opacity 
k of the gas and on the radiation flux <p coming from the core of the star as 



c 



(6) 



The ratio f/g of radiation pressure acceleration to surface gravitational acceleration 
determines how much mass loss takes place. The opacity k is an increasing function of 
metal content. In a population III star, opacity and thus mass-loss are expected to be 
mu ch lower than in a solar-metallicity star. 

Vink et al.l (1200 lb compute the mass-loss rate of stars for various metallicities and 
stellar masses. They find that a 6OM star with Z = 10~ 2 Z o loses about 10~ 7 M G /vr 
via radiatively driven winds. A Population III star living for » 10 7 yr contains even 
less metals, leading to a cumulative mass-loss of less than IM over its lifetime. The 
mass-loss is thus negligible and the mass stored in the star is fully available to produce 
an IMBH. 

Simulations show that non-rotating population III stars with masses in the range 
from 50 to 140 M Q collapse directly to a black hole, while b etween 140 and 260 Af n , 
explo sion as a Pair -Instability SuperNova (PISN) is expected (IHeger & Wooslevf2 002: 
Glover et al.ll2008l) . Masses higher than 260 M again form a black-hole. Formation 
of IMBHs from population III stars predicts the existence of a zone of avoidance in 
mass corresponding to the mass range where a PISN occurs, as this kind of supernova 
is expected to leave no remnant. 



2.2 IMBH formation from VMSs obtained by stellar merging 

The runaway merging of high-mass main-sequence stars in a den s e you ng cluster is 



a candidate mechanism for IMBH formation. Early on, IColgatd d 1967b studied the 
issue analytically, and in the last d ecades, extensive numerical simulations of young 
clusters have become widespread dFregeau et al.l20 06). These studies show that a Very 
Massive Star (VMS) forms if the lifetime of massive main-sequence stars (» 3 Myr) 
is longer than the typical time it takes for the runaway stellar merging to take place. 
Collisional stellar coalescence has to be fast enough to occur before the explosion of the 
first supernovae, which cause a sudden expansion of the cluster, reducing the number 
of physical collisions between stars. If a VMS is produced, its evolution likely leads to 
an IMBH remnant. Some high- mass stars that ha ve been observed in the local universe, 
such as the Pistol star (e. g., see iFigeret al.lll998b could be VMSs produced from such a 
stellar merging process (Freitag et al. 2006). In the following we will look more closely 
to the simple argument of lColgate ( 1967b . 
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2.2.1 A simple model of runaway merging 

Two stars of radii R\, R2 and of masses M\, M%, with typical relative velocity at infinity 
Voc, have a cross section for physical collisions in the form: 

,/ 2G(M l +M 2 )\ 

which differs from the purely geometrical cross-section because of the velocity-dependent 
expression enclosed in the rightmost bracket, which takes into account gravitational fo- 
cusing. Basically, stars bend each other's orbit via gravitational attraction, resulting in 
an enhanced collision cross-section. The ratio between the relative velocity at infinity 
and the escape velocity at contact 



2G(M l+ M 2 ) 

VE ~i c«i+*2) (8) 

determines whether gravitational focusing is important. Eq.|7]can be rewritten as 



cr = 7T(Rj + R 2 f 



V F 



2 \ 

(9) 



so the two limiting cases are Vg <K Voo, where gravitational focusing can be neglected, 
and Ve » Voo where gravitational focusing is dominant. 

For simplicity we will neglect any mass lost in stellar collisions, so all mass ends 
up in the star resulting from the collision. Let n be the number density of stars in the 
cluster and let us use the velocity dispersion (v) of cluster stars to estimate Voo- A given 
star collides with 

^ = ncr(v) (10) 
at 

other stars per unit time. If m is their average mass, the star's mass evolves as: 

dM 

= mncr{v) (11) 

at 

where M is the mass of the accreting star. 

The radius of a star scales approximately as R oc M 13 . Escape velocity at the star's 
surface scales as V§ = M/R oc M l ~P, so it increases with mass provided ft < I. For a 
massive star, gravitational focusing dominates and mass accretion is a positive feedback 
process. Under these assumptions, the cross section for a massive star becomes 

2 2GM 2nGM l+ P 

cr » jtR 2 5- « 12) 

RV 2 V 2 

where the mass (and radius) of smaller stars has been neglected. Using Eq. [TTJ we 
obtain: 

dM 2nGM l+ l 3 

— =mn (13) 

dt Voo 
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At time t — the (soon to be a) VMS had a mass Mo and a radius Rq, so: 



where 



so 



d_M _ Mo(M_Y +P 
dt t \M J 



{V) (15) 



2nGnmR{) 



which at time 



M(t) = M (l-^ ' (16) 

toiv = T ~ (17) 

formally diverges. A realistic value of /3 is 1 /2, so the divergence takes place in a finite 
time. For (3 — > the growth would instead become exponential. 
A VMS can thus be formed if 

tDiv < tMS ~ 3Myr (18) 

where tus is the typical lifetime of massive main-sequence stars, i.e. the time it takes 
for supernova explosions to take place. 

Another way to express the collisional merging timescale r is 

T = —T cr (19) 
«() 

where Rc is the radius of the young star cluster, T cr is its crossing time and 7?o is the 
typical radius of a star in the cluster. This naive timescale is actually too long for 
runaway merging to occur in a cluster with the typical crossing time and core radius 
observed today. To study the infall process and the subsequent runaway merging, thus 
going beyond the simple model described above, numerical approaches are needed, 
and they must necessarily take into account core-collapse and the very high densities it 
attains. 



2.3 IMBHs formed from merging of lower-mass black holes 



Miller & Hamilt on ( 2002) propose that dynamical interactions between binaries con- 
taining black-holes lead to repeated merging, which may result in the formation of a 
w 10 3 M Q IMBH if the process continues for billions of years, as may be the case in 
old star clusters, such as galactic Globular Clusters (GCs). Three-body interactions are 
crucial for tightening the binaries, bringing about coalescence, but they risk to push 
the system out of the cluster, thus stopping the merging process. Therefore the seed 



black hole must be massive enough (at least as 50 M Q according to Miller & Hamilton 
d2002h l to avoid expulsion. 



In a GC, the scale binding energy per unit mass is 



U> = J°jpL (20) 
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where a is a dimensionless shape factor. By setting 

Voo = ^lU b (21) 

and, taking the typical values of a quite concentrated GC, Mac — 10 6 M , Rcc = 1 P c , 
we obtain 

2GMcc =ax 92km/s (22) 
Rcc 

whereas for a looser cluster, M cc — 10 5 M e , Rcc = 4 pc, and 

Voo = a x 21km/ s (23) 

Thus, the typical escape velocity of a GC is about « 50 fcra/s. 

Let us consider a BH binary, where the primary has mass M and the secondary has 
mass m <K M. In an interaction with a third object of mass « m the binary hardens (i.e. 
increases its binding energy) by 

AU = ri^U (24) 

where 77 is typically « 0.2 dOuinlanll996h . A fraction » m/M of the eneigy At/ lost by 
the orbital motion is transferred to the center-of-mass motion of the binary because of 
momentum conservation, so its overall kinetic energy increases by 

Ar = ^t/ (25) 

To unbind the binary from the cluster, an increase of order the escape energy E esc is 
needed 

Ar = E esc = l -Mv 2 esc (26) 

so 

n W u = 5 M *" (27) 

and solving for U we obtain the minimum (in modulus) binding energy for which the 
binary is liable to be ejected due to a three-body interaction 

1 

2m 2 7] 



Uej = -r-^ (28) 



U e j must now be compared to the minimum energy U co for which the binary co- 
alesces, due to gravitati onal wave emission. The time it takes the binary to coalesce 
with this mechanism is dMiller & Hamiltonll2002h 

5c 5 a 4 , 2 .7/2 

c uw^m^ ' { y) 

for a binary orbit with semi-axis a, eccentricity e and reduced mass /u = mM/(m + M) » 
m. Equation |29] depends on orbital eccentricity, which varies in time, but we can take 
an average value (<?). So, using Eq.|29] we obtain 



_ GMm 



5GC d-(ef) V2 



128 



I/4 M 1/2 m 3/4 



(30) 
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and t c must be shorter than the time t e between two subsequent three-body encounters, 
otherwise the binary gets a new three-body kick before being able to merge. For t c = t e 
we obtain: 

5Gc 5 . , „ 7f2 \ l/A M^m 3 ^ 
(e) ) 



128 



■(1 



1/4 



Now, if U co < U e i, merging can take place before ejection. Thus 



5Gc 5 2 7/2 

l28- (1 - <e)) 



1/4 M" 2 m 3 ' 4 MV 



J/4 



2m 2 i] 



and by setting 



we obtain 



k = 2r] 



— (l-(e) 2 ) 
128 1 7 



7/2 



1/4 



M > k 2/5 



Gc m 



An alternative way to write this condition is 



M > k 2 ' 5 



(v eS c/c) 



1/10 



1/10 



(31) 



(32) 



(33) 



(34) 



(35) 



where r s is the Sc hwartzschild radius of the s tar of mass m, and the product v esc t e was 
indicated with A e . iMiller & Hamilto n (2002) adopt r\ » 0.2 and (e) « 0.7, obtaining 
0.1. If m = 10M o , v esc = 50 fcm/s and f e = 10 6 yr, Eq.[34lbecomes 



M > 70M G 



(36) 



This is approximately the threshold mass for IMBH growth through stellar-mass black 
hole merging. It is too large to come directly from stellar-mass black holes, so black- 
hole coalescence appears to need a seed black hole produced from some other process 
in order to work. On the other hand, once the accretion process is started, the IMBH 
grows quite fast, acquiring a 10 M through merging events every =s 10 yr, which 
over the lifetime of a GC (» 10 10 yr) can make it grow to « 10 5 M in mass, if enough 
stellar-mass black holes are present. 



3 Observational evidence of IMBH existence 

Evidence of the existence of IMBHs is still controversial. We will review two groups 
of phenomena: effects on the dynamics of stars in a cluster, and ultra-luminous X-ray 
emission. Several more phenomena are linked to IMBHs and could in principle yield 
a detection in the future: we sketch a discussion of these, but we do not add too much 
detail, as it would be premature. 
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3.1 Stellar dynamics in clusters 



The best way to detect an IMBH in a star cluster would be to observe the orbits of 
stars bound to it. This would allow to measure the black hole mass, as in the case 



of the super-massi ve black hole in the center of the Milky Way (see Gille ssen et al 
l2009bt iGhez et al.ll2008t billessen et ai1l2009at iGhez et alJl2005l Il998h . resulting in 
an uncertainty of m 10%, which is remarkably small by astrophysical standards. Sadly, 
no study of this kind has been completed with a positive result in any star cluster. This 
is understandable if we recall that the surroundings of the black hole in the galactic 
center were observed for more than 10 years, using state-of-the-art imaging technology 
(10 m telescopes with adaptive optics at KECK). An observational program of this 
kind on Galactic GCs, which are approximately 200, would literally require millennia. 
Moreover, cluster cores are extremely crowded (i.e. the images of several different stars 
will overlap) and the stars we need to observe are probably faint, main sequence stars, 
unless a red giant or an horizontal-branch star gets trapped in the IMBH potential well 
by chance. 

The methods used to detect IMBHs must then be indirect, at least at this stage. De- 
tection methods based on stellar dynamics can be divided between those which probe 
the zone of the cluster directly affected by the presence of the IMBH, known as the 
IMBH influence sphere, and those who do not need to. The former are more affected 
by crowding in the cluster core and by small number statistics, because the influence 
sphere usually contains few stars. 



3.1.1 The influence sphere 

Let the IMBH mass be Mbh and the cluster mass be Mqq- The size of the region where 
the IMBH strongly affects the dynamics of stars may be calculated by equating the 
escape velocity v e at a distance Rinf from the IMBH to the velocity dispersion <x in the 

cluster core: 

2GM BH 

~R = °" (3?) 

t<INF 

and 



GM C 



Rgc 

where Rqc is a scale radius for the cluster. So, neglecting factors 

Rinf Mbu 



Rgc Mr, 



(38) 



(39) 



for a typical globular star cluster, Rgc ~ 1 pc. If we take Mbh/Mcc - 10~ 2 , then 
Rinf = 10~ 2 pc. The distance of the nearest globular star clusters from the Sun is on 
the order of several kpc At a distance of * 1 kpc the angular size of the influence 
sphere is a 2 arcsec For nearby clusters it is thus feasible to probe the influence 
sphere with present-day telescopes (e.g. Hubble Space Telescope). The volume of the 
influence sphere is 



3 

Mcc, 
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The number density of stars in the cluster core is a factor of 77 « 10 -f 100 larger than 
the average density, so the fraction of stars within the influence sphere is 



NlNF I M BH 

= n\ 



GC 



M, 



GC 



and given that Ncc ~ 10 5 - 10 6 for a typical globular star cluster we obtain 



N, 



INF 



10- 10 2 



(41) 



(42) 



which is a small number if we hope to obtain a smooth density profile. This is one of 
the major hindrances encountered in revealing IMBHs using the dynamics of stars in 
the influence sphere. Still, some interesting results can be obtained, as shown in the 
following section. 



3.1.2 Evidence within the influence sphere 

Peebles! (1 19721) and lBahcall & Wolfl d 19761) were amongst the first to predict the forma- 
tion of density cusps in star clusters containing a central massive object. Both find a 
phase space distribution function of stars near the central mass in the form 



f{r,v)=K(-EY 

where E is the energy per unit mass (assuming that all stars have the same stars) 

v 2 GM BH 

E 

2 r 

and r is the distance from the central object. The associated density profile is 



n(r) 



where 



I p — I (1 - x) x 2 dx 
Jo 



(43) 



(44) 



(45) 



(46) 



so a pow er-law cusp in the three-dimensional density profile is predicted, with exponent 
p + 3/2. Bah call & Wolfl (119761) finds p — 1/4 based on an analytic argume nt. This 
result is confirmed by simulations based on different numerical approaches dShapirol 
1985|:lFreitag & Benzl2002l:lBaumgardt et al.l2004lPreto et al.l2004tlAmaro-Seoane et al 



20041) . Observing a cusp in the core of a GC may thus be an indication of the presence 
of an IMBH, but some caveats apply. Cusps may form due to other physic al processes 



such as core collapse in a cluster containing stars with different masses ( Tre nti et al 



2010h . On the other hand, a cusp may easily go undetected if the center of the cluster 
is incorrectly determined or if the resolution of the o bser vations is limited. Notwith- 
standing these limitations Novola & Gebhardtl (12006b and iNovola & Gebhardtl (12007b 
were able to reveal central luminosity density cusps in a fraction of the sample of GCs 
they observed using Hubble Space Telescope. 



13 



A better grip on the central dyna mics of cluste r s can be obtai ned if spect r oscopic 
(i.e. kinematic) data are also used. iNovola et al. (2008a) and INovola et alj (2008b) 



obtained kinematic information for the inner 5 arcsec of a> Centauri using the Gemini 
telescopes, resulting in a claim to the detection of a 4 x 10 4 M black hole in its center. 
The point here is that the observed stellar velocities increase near the center (i.e. a 
velocity dispersion cusp is observed). A central rise in the velocity dispersion argues 
for the presence of a dark mass which deepens the cluster potential well in the cen- 
ter. However, the kinematic information form spectroscopy is limited to line-of-sight 
velocities, which might not tell the whole story. A strong anisotropy in the central 
pressure tensor could account for the observed rise, as line-of-sight velocities could be 
much higher than velocities on the plane of the sky. Moreover, centering is still a del- 
icate issue here, and the mechanisms which could generate a density cusp (see above) 
are likey to generate a velocity disp ersion cusp as well 



At variance with the claim of INovola et a l. (2008a), Ivan der Marel & Anderson 



(2010) Ind an upper limit of 1.2 x 10 4 M to the mass of an IMBH in lo Centauri. 



A definitive solution to the puzzle has not been found yet. 
3.1.3 Evidence outside the influence sphere 

Evidence depending on resolving the influence sphere is limited by the low number 
of stars contained in it. Methods which rely on the global cluster structure to test for 
IMBH presence are free from to these shortcomings. The density profile of a GC, either 
three-dimensional or sky projected, is characterized by two characteristic radii and the 
overall luminosity. The central part of the cluster is almost constant in density, up to 
a core radius, R c , from which the density starts to decline. So R c measures the size of 
the core of the cluster. The overall size of the cluster, instead, is characterized by R/,, 
the radius conta i ning h alf of the cluster mass. 



Trenti et al.l d2007l) use N-body numerical simulations to find that, regardless of 
the initial conditions, clusters with an IMBH of about 1 % the total cluster mass and a 
realistic content of binary stars (about 10%) settle to a ratio of the (three-dimensional) 
core radius to the (three-dimensional) half-light radius of about 

^ * 0.3 (47) 

after becoming sufficiently dynamically old. This is a; 3 times as large as the ratio 
obtained for clusters without an IMBH. The core of a cluster containing an IMBH is 
larger in proportion to the overall size of the cluster. 

Let us provide a justification of the above result from the physical point of view. We 
mentioned in the previous sections that a GC loses mass and energy as stars evaporate 
from it. Evaporation is due to two-body relaxation, i.e. stars exchanging momentum 
and energy in gravitational encounters. Two body relaxation promotes thermaliza- 
tion, i.e. the cluster phase space distribution function is driven towards a maxwellian 
over timescales of about one to tens of Gyr, depending on the cluster radius and num- 
ber of stars. A maxwellian, however, has a tail with arbitrarily high velocities. The 
stars ending up in such a tail must be ejected from the cluster, removing energy and 



14 



mass from it. The evolution in time of the half-mass radius is influenced by this evap- 
oration process as follows ( the following ana lytical argument is loosely inspired by 



Vesperini & Cherno ff ( 1994) and lHenonl d 19651) ). 



The potential energy of the cluster is 



1 v~i Gmitn j 

2 2^ (48) 



U 2^ , 

where m, is the mass of the z'-th star, and r, ; is the distance between the z'-th and the y'-th 
star. The virial theorem holds: 

2T + U = (49) 

so 

E = T + U = | (50) 
We introduce an effective radius of the cluster such that 

R eff 

where M = Yii m i- Note that in general R e ff + R c + Rh. So 

2M 2 _ 2N 2 

R ?ff ~ y mtnj ~ — T ( 52 ) 

where is the total number of stars. The last equality holds approximately for equal- 
mass stars. The radius R e jf can be very different from Rh, as if any — > 0, the 
denominator diverges and R e ff — > 0. In other words, if stars generate subclusters 
(e.g. binary or triple systems), their binding energy becomes extremely high, giving a 
significant contribution to the total binding energy of the cluster. The ratio 

n = ^ (53) 

measures the amount of subclustering present in the system: the larger 77, the less 
subclustering is present. From Eq.lBTIand Eq.|50]we obtain 

E -^- k -^L (54) 
E M R eff 

where the dot indicates a derivative with respect to time. Eq. l53limplies 

^IL = i + ^ (55) 

R eff n R h 

SO 

KJL-_ 2 ™A^ (56) 
R h M E 77 
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The cluster loses mass due to evaporation, so M/M < 0. It loses energy as well, but its 
binding energy is negative, so E/E > 0. Then 2M/M — E/E < 0. If no substructure 
forms, t) = and the only way the cluster has to satisfy the virial theorem and conserve 
energy is to globally shrink, reducing R/,. Shrinking increases the binding energy, com- 
pensating for the energy lost to evaporation. If instead 77/7/ is negative, i.e. substructure 
forms, the cluster can keep energy balance without shrinking. Lower values of 77 cor- 
respond to the formation of hierachical subclusters, so binary stars form (or harden) 
in order to increase the cluster binding energy while allowing for a constant or even 
expanding R/,. This makes physical sense, as binaries can form only via three-body in- 
teractions (the third star acts as an energy sink). Three body interactions need an high 
density of stars to happen in significant numbers. An initial shrinkage of the cluster 
makes the central density rise, preparing the conditions for three-body interactions to 
take place and for binaries to be generated, or hardened if they are already present. 

In simulations, the IMBH is usually found in a binary with a stellar-mass black- 
hole, so it acts as a source of energy, giving a significant contribution to the substructure 
term 77/77. Sources of energy in the form of substructure (several binaries or an high- 
mass binary with the IMBH as principal) can halt the collapse of the cluster and sustain 
a slight expansion in Rf, as well as a large core with high R c /Ri,. This explains the result 
expressed by Eq. [47] 

Observationally, R c /Ri, can be used as a very preliminary method to select GCs 
which are eligible as IMBH hosts, but several difficulties are present. Sky-projected 
quantities only are measurable in an observed GC. The half-mass radius is not directly 
observable, as only luminous matter can be used to compute it, and the dark remnants 
(neutron stars, stellar mass black holes, white dwarfs) are typically more massive than 
visible stars and tend to segregate to the center of GCs, thus leading to an overestimate 
of Rh as in general the light profile does not trace the mass profile of the cluster. This is 
a major hindrance, as the relative amount of stellar-mass black holes and neutron stars 
with respect to the visible stellar component is poorly constrained. 



3.2 Ultra Luminous X-ray sources (ULXs) 
3.2.1 Eddington limit 

A source of radiation of mass M held together by its own gravity cannot have an arbi- 
trarily high luminosity. A particle of mass m and charge q at distance r from the source 
feels gravitational attraction to the source 

77 GMm «T» 
F g = —j- (57) 

and repulsion from it due to radiation pressure 

F >-i^ (58) 

where <x = <j(q) is the cross section of the particle to the radiation from the source, and 
L is the source luminosity. Equation|58]holds only in the case of spherically symmetric 
emission. 
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By requiring that F r < F g so the particle is globally attracted to the source we 
obtain 

GMmAnc 

L < = L EDD (59) 

so for any mass M of the source there exists a critical luminosity Ledd at which radia- 
tion pressure would destroy the source. We can take <jj e (the Thompson cross section) 
as the electron cross section 

8tt 2 %nq\ 

where r e is the classical electron radius, q e and m e the electron's charge and mass. 
Instead for a proton we obtain 

aTp = 3m^- (61) 

Thus a-jp <k cr Te , because m e <k m p . As an estimate for hydrogen, we use m p 
1 .67 x 10~ 24 g for m and crj e — 6.65 x 1CT 25 cm 1 for cr, obtaining 

Ledd = — X 1.26 X 10 3 Vg/s (62) 

Mq 

as the maximum allowed luminosity for the source not to be destroyed by its own 
radiation pressure. 

3.2.2 Observing and explaining ULXs 

An unresolved X-ray source (i.e. a source that appears point-like in observations) is 
defined as ultra-luminous if 

L x > I0 39 erg/s (63) 

i.e. if its luminosity in the X-ray band exceeds the Eddington luminosity of an a 1OM 
mass. The centers of galaxies often display X-ray activity, but they are usually excluded 
from this class of objects. Thus only those sources which are off-center with respect to 
the parent galaxy are counted among the ULXs. 

High-luminosity X-ray sources began to be discove red in galaxies thanks to the 
Einstein Observatory mission in 1978 (Fabbianq[l989|). ROSAT ASCA, Chandra 



and XMM-Newton followed, resulting in a wealth of data on the X-ray sky. ULXs are 
widespread in nearby galaxies, but a satisfa ctory understan ding of the underlying phys- 
ical mechanisms has not been attained yet ( Fabbia nol2006h . To achieve luminosities as 
high as 10 39 erg/s, the only plausible physical mechanism proposed so far is accretion 
of matter onto a compact object. However, by definition, ULXs are impossible to ex- 
plain by isotropic, sub-Eddington emission from an accreting stellar mass black hole 
of M < 1OM . 

The only remaining viable explanations are that either ULXs are beamed sources 
or somehow achieve super-Eddington emission, or a new class of accreting objects, 
that is IMBHs, is the s olution to this observational riddle. We point the reader to 



Colbert & M iller (2005) for a review of the observational evidence supporting the latter 
point of view. In the following we briefly review the most important piece of evidence, 
that is spectroscopic data. 
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3.2.3 Spectroscopy and the MCD model 

The Multi Color Disk (MCD) model was first used to fit the X-ray spectra of Galactic 
sources emitting at lower luminosities than 10 39 erg/s. It stems essentially from the 
superposition of a blackbody spectrum and the so-called multi color disk spectrum 
(Mitsu da et alJll984t iTakano et al.l [l994). Galactic X-ray sources are neutron stars or 
stellar-mass black holes accreting mass from a companion star. The accreting matter 
forms a disk orbiting the neutron star or black hole and increases its temperature as it 
nears its surface. The blackbody component corresponds to the surface of the accreting 
object, and the multi color disk component is a sum of blackbody sources at different 
temperatures intended to model the accretion disk. 

A particle of mass m in a newtonian orbit of semi-axis r around the accretor mass 
M has an energy 

GMm 

E = -— — (64) 
2r 

so, a change of its semi-axis by an amount dr influences the energy as 

ip GMm 

dE = — dr (65) 

2r z 

Such a change can be brought about by a perturbing force, most likely friction due to 
gas viscosity in our case. 

The energy Ed deposited per unit time by the accreting matter in the region between 
r and r + dr scales as Ed oc dr/r 2 . The energy Ej radiated per unit time, on the other 
hand, scales with the surface of the emitting region, i.e. is proportional to rdrT 4 , so 

-% oc rdrT 4 (66) 

which, integrating, yields 

Tocr 3 ' 4 (67) 



so for a face-on disk at distance D from the observer we obtain a spectrum 
~ r » 8 7rr? r T '(TV n/3 _ 

D- ' 

where 



1 r r ° %-nr 1 , rT,i T \-^l^ 



(68) 



E 3 dE 

I{E,T)dE = A {emT _ i) (69) 

is the Planck distribution, rj is the innermost radius of the disk, nearest to the accreting 
object, ro the outermost radius and Tj and To the corresponding temperatures at these 
radii. Notice that ro s> rj and Tj > To- 

The parameter 77, which is thus obtained by fit ting the spectrum, depends on the 
mass of the accreting object (Makishim a et al h00Ch : 



Tj oc AT 1/4 (70) 



so 



lower temperatures of the inner edge, correspond to higher masses of the accreting 
iect. 



object. 
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So high-luminosity ULXs with spectra that are well fit by a MCD model with low 
77 are good IMBH candidates. Typical temperatures obtained from Galactic X-ray 
sources associated with stellar-mass accreting objects are in the 0.4 - 1 .0 keV range, so 
the expecte d temperature for an obje ct « 100 times more massive will range from 0.1 
to 0.3 keV. IColbert & Milled d2005l) (Table 1) report several ULXs having 7> in this 



range and a luminosity in excess of 10 erg/s, which are potential IMBH candidates. 
3.2.4 ULXs and star-forming regions 

As a last remark on ULXs, we notice that they appear to be associated to active star- 



forming regions (Humphrey et al. 2003; Swartz et al. 2004; Liu & Bregman 2005), in 



partic ular interacting galaxy pairs: e.g. the Antennae (N GC 4038/4039) dFabbiano et al 



200 lb and the interacting galaxy pair NGC 4485/4490 dRoberts et alj|2002h . If ULXs 



are really associated to IMBHs this feature would support the runaway merging model 
of IMBH origin, as young clusters with high-mass main sequence stars are found vir- 
tually only in star-forming regions. 

3.3 Other phenomena possibly related to IMBHs 

Other phenomena have been regarded as possible signs of the presence of an IMBH 
in a star cluster. Some of them are reviewed in the following, with the caveat that the 
listing may be incomplete due to the large amount of work published on the subject. 

3.3.1 Millisecond pulsar timing 



D'Amico et alj ([2002) use the timing of three Milli-Second Pulsars (MSPs) in the GC 
NGC 6752 to obtain lower limits to their line-of-sight acceleration due to the cluster 
potential. This implies a lower limit of about 10 M /L for the central mass-to-light 
ratio, equivalent to 1.3 x 10 4 M in dark objects within the innermost 0.15 pc around 
the cluster center. IMBH of this mass would be a solution, but such a central increase 
in the mass-to-light ratio can be explained away by dark remnants (stellar-mass black 
holes and neutron stars) which undergo mass-segregation and accumulate in the center 
of the cluster potential well. 

3.3.2 High-velocity stars 



Trenti et alj d2007l) shows, among other results, that a population of hard binaries inter- 



acting with an IMBH can produce several stars which escape the cluster wi th a velocity 



of hundreds of km/s, much higher than typical cluster escape velocities. Hillsl ( 1988 ) 



proposed to use high-velocity stars (possibly) eject ed from the galactic c enter as a 
test for the presence of a super-massive black hole. EE & Tremain e (2003) estimates 



the number and the velocity of high-velocity stars produced by various mechanisms 
(gravitational encounters of single stars, tidal breakup of a binary star, ejection by a 
binary black hole) analytically in quantitative detail. In principle these estimates can 
be adapted to star clusters, yielding a potential method for detecting IMBHs. 
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